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Abstract 

Recently, Galley [Phys. Rev. Lett. 110, 174301 (2013)] proposed an initial value problem formulation of Hamilton’s 
principle applied to non-conservative systems. Here, we explore this formulation for complex partial differential 
equations of the nonlinear Schrodinger (NLS) type, examining the dynamics of the coherent solitary wave structures 
of such models by means of a non-conservative variational approximation (NGVA). We compare the formalism of the 
NCVA to two other variational techniques used in dissipative systems; namely, the perturbed variational approximation 
and a generalization of the so-called Kantorovich method. All three variational techniques produce equivalent equations 
of motion for the perturbed NLS models studied herein. We showcase the relevance of the NGVA method by exploring 
test case examples within the NLS setting including combinations of linear and density dependent loss and gain. 
We also present an example applied to exciton polariton condensates that intrinsically feature loss and a spatially 
dependent gain term. 
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1. Introduction 

Variational methods are commonly used to describe 
the dynamics of nonlinear waves in nonlinear optics, and 
atomic physics i i a i- These methods rely on a 
well-informed ansatz substituted in the Lagrangian or 
Hamiltonian formulation of a complex, infinite dimen¬ 
sional system. This ansatz reduces an original partial 
differential equation (PDE) model to a few degrees of 
freedom establishing equations describing the approximate 
dynamics appropriately projected into the solution space 
spanned by the ansatz. The variational approximation 
(VA) method projects the high-dimensional dynamics to a 
low-dimensional dynamical system for the time-dependent 
parameters that encapsulate the qualitative and quantita¬ 
tive behavior of the original complex system. It is im¬ 
portant to note that there is an intrinsic drawback of 
VA methods in that they have a strong restriction when 
projecting the infinite-dimensional dynamics of the orig¬ 
inal PDE to a small hnite-dimensional ansatz subspace. 
This OTojection is known to potentially lead to invalid re¬ 
sults [6(1 , a feature which is naturally expected (given the 
large reduction in the number of degrees of freedom) when 
the full PDE dynamics ceases to be well-described by the 
selected ansatz. Nonetheless, there have been some efforts 
to control the corrections of the VA to increase the ac¬ 
curacy of the results Q- Fundamentally, the variational 


method relies on the existence of a Lagrangian or Hamil¬ 
tonian structure from which the Euler-Lagrange equations 
can be derived. This prerequisite limits the application of 
the variational approach to conservative, closed systems. 

The recent work by Galley [H offers a new perspective 
to the classical mechanical formulations by recognizing 
that the Hamilton-Lagrangian formulation has the key fea¬ 
ture of being a boundary value problem in time although 
it is used to derive equations of motion that are solved 
with initial data. By treating the extremization problem 
as an initial value problem, variational calculus can be ap¬ 
plied to non-conservative systems. Although Galley’s pro¬ 
posal was for classical mechanical, few-degree-of-freedom 
systems, i.e., systems described by ordinary differential 
equations (DDEs), it paved the way for its application to 
dispersive nonlinear PDEs. The relevance and validity of 
the extension of Galley’s method to nonlinear PDEs was 
showcased in the recent work of Ref. Q where examples 
cases based on a 'PT-symmetric sine-Gordon and mod¬ 
els were presented. This work recognized that not only 
can the method be utilized for infinite degree-of-freedom 
systems, but its Lagrangian underpinning enables a varia¬ 
tional approach to be developed for non-conservative PDE 
systems. Ref. @ illustrated very good agreement between 
the original PDE dynamics and the reduced ODE model, 
obtained from the corresponding non-conservation varia- 
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tional approximation (NCVA). 

In this manuscript we present the extension of this 
NCVA method for the specific, but broadly impor¬ 
tant/applicable case of the nonlinear Schrodinger (NLS) 
equation. The latter model and its variants are of princi¬ 
pal interest to applications from optical physics @ , atomic 
physics 0 and other areas of mathematical physics [lll |. 
not only in their conservative, but also in dissipative vari¬ 
ants of the model |l2|. The details of our presentation 


are as follows. Section [2] presents the variational formal¬ 
ism. For completeness, the section starts with a short 
review on Galley’s method as originally introduced for 
finite-dimensional Hamiltonian systems and extended to 
complex non-conservative forces. Then, we extend this 
methodology to the NLS equation and compare it to the 
perturbative variational approximations, examined previ¬ 
ously. We, in fact, prove that these different methods, 
within the NLS setting, are equivalent. In Sec. [3] we 
present a few examples of NLS settings including non- 
Hamiltonian terms. We use a combination of linear and 
density dependent loss and gain and showcase an exam¬ 
ple in the real of polariton condensates that are out- 
of-equilibrium as they intrinsically contain loss and gain 
terms. Finally, in Sec. 0] we present our conclusions and 
discuss a few avenues for further exploration. 


2. Non-conservative 
Formalism 


Variational Approximation 


2.1. Non-conservative Approximation for Finite- 
Dimensional Systems 

Let us discuss the non-conservative variational approach 
(NCVA), based on the non-conservative variational formu¬ 
lation introduced by Galley in Ref. [i[. Our aim here will 
be generalize the latter to include complex variables and 
complex non-conservative forces. Subsequently, we will 
use the NCVA to study the dynamical characteristics of 
nonlinear waves in the perturbed NLS model. In Ref. [l|. 
Galley highlights that the time-symmetric and conserva¬ 
tive dynamics in Hamiltonian systems is due to the bound¬ 
ary value form of the action extremization problem. As 
such, the original formalism based on extremization of an 
action functional is predicated on the conservative nature 
of the system. Therefore, a direct application of action 
extremization when the system includes non-conservative 
terms is, by construction, not viable. As Galley points 
out, in simple classical mechanics cases where the dissipa¬ 
tion forces are local in time and linear in the velocities, it 


is possible to employ Rayleigh’s dissipation function [13 1. 
However, this formulation cannot be applied to systems 
with more general dissipative forces such as nonlocal or 
nonlinear ones. 

Therefore, in order to treat systems with general forces. 
Galley proposed to consider the extremization problem as 
an initial value problem instead. The method is based on 


considering two sets of variables, qi and q 2 , and applying 
variational calculus for the non-conservative system pro¬ 
vided q 2 = qi after the variation. Gonsider the path q{t) 
passing through the (fixed) boundary values (f at t = ti 
(initial) and qj at t = tf (final). Suppose that the sys¬ 
tem trajectories are described by the following set of N 
generalized coordinates and velocities: q = {qi}^i and 
q = {qi}fLi- Let us now double both sets of quantities, 
q ( 01 , 92 ) and q ( 91 , 92)5 and parametrize both coor¬ 
dinate paths: 


9i.2(t,e) = 9i.2(i,0) -he 771 , 2 ( 0 , (1) 


where 9*1,2 (t, 0 ) are the coordinates of two stationary paths 
(e <C 1 ) and f]i, 2 {t) are arbitrary virtual displacements. 
The following equality conditions are required for varying 
the action such that the endpoints coincide: 771 , 2 (^ 1 ) = 0 , 
9i(t/,e) = 92(t/,e), and qi{tf,e) = q 2 {tf,e). Therefore, 
the equality condition does not fix either value at the final 
time. After all variations are performed, both paths are 
set equal and identified with the physical path, q{t), the 
so-called physical limit (PL). 

The total action functional for 91 and 92 is dehned as 
the total line integral of the Lagrangian along both paths 
plus the line integral of a function TZ which describes the 
generalized non-conservative forces and depends on both 
paths {qu}l^i- 


= / dt£(9i,9i)- / dt£(92,92) 


-j- / dtlZi^qrfi j qrfit') ^ 


dt 


-^(91,91) --^(92,92)+'^( 9 ™, 9 ™,i) ( 2 ) 


The above action defines a new Lagrangian: 

Hqn,qn) = £( 91 , 9 * 1 ) -£( 92 , 9 * 2 ) -l-7^(9„,9„,t). (3) 


If TZ is written as the difference of two potentials V{qi) — 
H( 92 ), i.e., under the presence of conservative forces, then 
it is absorbed into the difference of two (suitably redefined) 
Lagrangians, leaving TZ effectively to be zero. On the 
other hand, a nonzero (nontrivial) TZ corresponds to non¬ 
conservative forces and couples the two paths together. 

For convenience, following Ref. [i|, we make a change 
of variables to 9 + = (gj -1- 9*2)72 and 9 _ = 91 — 92 - In 
this way, 9 _ —^ 0 and 9 + —^ in the physical limit. The 
conjugate momenta are then q± = dA/dcffi and the paths 
are parametrized by q±{t,e) = q±(t,0) -\- erf±(t). There¬ 
fore, the new action is stationary under these variations if 
\dS[q±]/de\e=o = 0 for all 77 ^. 

In the original 91,2 coordinates the equations of motion 
may be written, using the Euler-Lagrange equations, as 
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dTTi^ 2 /dt = dA/dqi ^2 with 7ri_2 = (—l)^’^9A/c}(fi^2- Simi¬ 
larly, in the q±, the Euler-Lagrange equations yield to the 
equations of motion. However, in the physical limit, only 
the dA/dq- = d7r-)-/c?t equation survives, such that the 
trajectory is defined by 


TT(q, ^ = 


dt 

with conjugate momenta 
Tr{q, ^ = 


’ dK ' 

dC 

'dTZ.' 

dq- 

PL 

dq- 


PL 


■ dA' 

_ dC ^ 

'dTZ' 

.dq-. 

PL ^9 

.dq-_ 


( 4 ) 


( 5 ) 


PL 


In the presence of conservative forces (i.e., 7^ = 0), the 
usual Euler-Lagrange equations are recovered. However, 
for non-conservative forces, a nonzero TZ modifies the 
trajectories as per Eqs. o and ©■ Since in our con¬ 
siderations below we are concerned with complex non¬ 
conservative forces, the non-conservative term TZ will be 
complex as well. The terms in TZ are precisely responsi¬ 
ble for coupling both the qi and q 2 paths to each other. 
It is important to note that the complex conjugate of 
the functional terms in C are necessary for solving the 
Euler-Lagrange equations when considering coordinates 
containing complex terms —which is precisely the case 
for the NLS under consideration. In the physical limit, 
only the Euler-Lagrange equation for the -I- variables sur¬ 
vives. Therefore, expanding the action in powers of q- the 
equations of motion follow the variational principle: 


Sq-{t) 


( 6 ) 


This also has the consequence that only terms in the new 
action that are perturbatively linear in (f_ contribute to 
physical forces. 


2.2. Non-Conservative Variational Formulation for Non¬ 
linear Sehrodinger Equation 

Let us now extend the NCVA formalism for the NLS 
equation. It is worth mentioning at this stage that the 
NLS equation is, arguably, a prototypical nonlinear PDE 
that supports travelling envelope waves. Namely, the NLS 
equation is the normal form for envelope waves 1^. As 
such, the NLS equation is a paradigm for wave propaga¬ 
tion and a universal model describing the evolution of com¬ 
plex field envelopes in nonlinear dispersive media 00. 
The NLS equation models a wide range of physical phe¬ 
nomena including hydrodynamic waves on deep water , 
nonlinear optical systems |^. [^ 21, [sll. heat pulses 
in solids 221, plasma waves U 2J, 251, and matter waves 
IB) 271; while it has also attracted much interest mathe¬ 


matically [13, [m, m, [ii [IB, m I. 


The one-dimensional (ID) NLS equation can be cast in 
non-dimensional form as |32| 

fwt + +5|upM = 0, (7) 

where u{x,t) is the complex field and g is the nonlinear¬ 
ity coefficient [g = -|-1 {g = —1) corresponds to an at¬ 
tractive (repulsive) or focusing (defocusing) nonlinearity]. 
This form of NLS corresponds to a Hamiltonian PDE. The 
Lagrangian density for this conservative NLS is [3B.l33.[33| 

£ = - uu;) -f (8) 

where (•)* denotes complex conjugation. Eor consistency 
of notation we will use calligraphic symbols (cf. C) to 
denote densities while their effective (integrated over all 
x) quantities we will use standard symbols. Namely 
L = C dx. In the presence of non-conservative terms 
ifP) that may depend of the field u, its derivatives, and/or 
its complex conjugate, the NLS takes the general form: 

iut + ]^u^j,+g\u\^u = F’. (9) 

Following Ref. [l[, two coordinates are introduced: ui and 
U 2 . In analogy with the finite-dimensional case of the pre¬ 
vious section, we construct the corresponding total La¬ 
grangian for these two coordinates: 

£t = — £2 + (10) 


where Ci = £(ui, .., t), for i = 1,2, represents 

the corresponding conservative Lagrangian densities for 
Ml and U 2 as defined in Eq. ([8|) and TZ contains all non¬ 
conservative terms that originate from the term T’ in 
Eq. ([7]). Therefore, by construction, the non-conservative 
part of the Lagrangian (fT0)l must be related to the pertur¬ 
bation term T’ in the NLS ([9]) by 


V = 


dTZ 


du* 


-J PL 


( 11 ) 


such that TZ = T’ u*_ -\- const, where the constant of inte¬ 
gration is with respect to u*_. As before, for convenience, 
M+ = {ui -\- U 2)/2 and u- = ui — U 2 are defined in such 
a way that at the PL u+ —>■ u and u- —^ 0. The cor¬ 
responding conjugate momenta are defined as in Sec. 12.11 
and thus, the equation of motion reads 


r STZ 

dtSu^ Su* _6u*L 


( 12 ) 


where S denotes Frechet derivatives. 

Through this method we recover the Euler-Lagrange 
equation for the conservative terms and all non¬ 
conservative terms are lumped into [5TZ/5u*_]^i,. It is cru¬ 
cial to construct an TZ such that its derivative with respect 
to the difference variable u*_ = u\ — U 2 at the physical limit 
gives back the non-conservative or generalized forces [P in 
Eq. (|ni)]. A similar variational formulation can be applied 
to other PDE (or ODE) models of interest. 
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2.3. Perturbed Variational Approach Formalism and 
Equivalence Proof 

Let us now compare the above methodology using 
the NCVA to the standard perturbed variational ap¬ 
proach 33 . Let us consider the non-conservative modified 
NLS equation Eq. ([H|) with the non-conservative general¬ 
ized force V = eQ, where e is a formal, small, perturbation 
parameter (|e| <C 1). In this manner, when e = 0, one re¬ 
covers the conservative Lagrangian ([8|) with corresponding 
Euler-Lagrange equations: 


^ ^ 

dp dt dp 


(13) 


(i.e., an entry of p) and remind the reader that equations 
with the symbol p denote a set of couple equations for 
each of the entries p in p. As in Galley’s prescription, for 
simplicity when evaluating in the physical limit, define the 
± coordinates for each parameter: = (pi +p 2)/2 and 

P- = (pi — P2) such that p_ = p*_, then we can show the 
NCVA is equivalent to the perturbed variational approxi¬ 
mation: 


P 


P + OO POO 

r dv 1 

/ Pdx = / 

—00 —00 

du*_ 


dx, 


(19) 


projected into the ansatz, such that 


where L{p) = f Cdx and £ = L[u{xp,p)\ is the conser¬ 
vative Lagrangian (|8]) evaluated on the chosen variational 
ansatz u bearing a vector of variational parameters p, on 
which the effective Lagrangian £ depends. For consistency 
we will now use a bar over quantities that are evaluated at 
the variational ansatz. To adjust for the presence of small 
non-conservative terms one needs to find the remainder of 

OLt d OLt - . 

dp ~ dt dp’ ^ ’ 


which is nonzero, where the Lagrangian Lt = L + 
conservative terms L and the non-conservative terms 
ApplyiM the (linear) perturbed variational approxi¬ 
mation [ 3 ^ Q yields, after obtaining the effective La¬ 
grangian from the Lagrangian density, the following per¬ 
turbed Euler-Lagrange equation: 


d dl dl 
dt dp dp 



^du* 

Q + 

dp dp 


dx. 


(15) 


The right hand side above is also equivalent to the result 
of the modified Kantorovich method [s^ yielding in the 
right hand side of Eq. m 



-^du -du* 

p* -- 

dp dp 


dx = 2Re 



T^-^dx. 

dp 


(16) 


Let us now establish the equivalence between the above 
perturbed variational formulation and the NCVA method. 
Given the non-conservative NLS m, where V is assumed 
complex, let us construct the functional TZ evaluated at 
the variational ansatz such that 


P = 


-00 [dp*- 




L 


“ ^-du* 
V ■ 


dp*_ 

dP 


dx, 


— 00 . 
00 




dp*_ 

du 


' dp*_ 
■du* 


-I PL 

- du- dp* 
V*-^+u- 


dp*_ 


' dp*_ 


dx, 


J PL 


dp* 


I 


( 20 ) 


since [mL]pl = [h-]pL = 0- The non-conservative inte¬ 
gral in the Euler-Lagrange equation derived in Eq. (1201) is 
equivalent to the perturbed variational approximation in 
Eq. (US]). 


3. Applications of the NCVA to the NLS 

3.1. Linear loss 

As a first example for the application of the NCVA, we 
use the focusing {g = -|-1) NLS equation with a linear loss 
term of strength e: 


iut + -Uxx + \u\'^u =-ieu. (21) 

In the absence of the linear loss (e = 0), the NLS (l?T]) ad¬ 
mits the well-known, bright, solitary wave solutions 
[^ . Therefore, in order to follow the effects of the lin¬ 
ear loss on the soliton, we choose a bright soliton ansatz 
with arbitrary height a, inverse width w, center position 
£ speed c, chirp b, and phase (j) as follows: 

UAix,t;p) = asech{w {x — -i-c(x-5)-i-<»^ 1 - 22 ) 


P = 


dn 


du*_ 


PL 


(17) 


Let us now require that for the evolution of the ansatz the 
variational parameters are real. Thus, the solution which 
satisfies Eq. dia and ensures real values for the parameters 
is: 

P = P{u±,u*,^,u±u,...) uL + C.C., (18) 

where c.c. stands for complex conjugate. For ease of no¬ 
tation let us denote by p a single variational parameter 


where the vector of time-dependent parameters corre¬ 
sponds to: p = {a,w,^,c,b,4>). For the conservative 
variational approximation, one can dehne an effective La¬ 
grangian L, with the expected Euler-Lagrange equations 
of motion. 

In the NCVA framework, the ui and U 2 ansatze are 
dehned as in Eq. (|2^ : 


Ml = UA{x,t;pi), (23) 

U2 = UA{x,t;p2), (24) 
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where each solution has its corresponding parameters: 
Pi = {ai,wi,^i,ci,bi,(t)i) and p2 = {a2,W2,^2,C2,b2,4>2)- 
According to the non-conservative variational method the 
Lagrangian is Ct = Li — C 2 + Ti- where 


of motion: 


/ 

a = —ab, 
b = 

TT^ 


— 26 ^, 


^1 = ^ -Mlwi.t) + (25) 

L 2 = ^{u2U2^t-U*U2,t) +^\u2,x\^ -^\U2\'^, (26) 

n = ie(u 2 ul-uiu;). (27) 

We then write the effective Lagrangian L = Lxdx = 
Lidx — L 2 dx + TZdx, for which Li and £2 
recover the same equations of motion as the conservative 
variational approximation. This yields the modified Euler- 
Lagrange equations: 


dp 


d 

dt 


dp 


dn 


dp- 


dx = 0. 


(28) 


PL 


After integration and simplification, the effective La¬ 
grangian is given by: 


L 


1 2 , 

—aiWi + 2' 
1 

+ 77 - 


WX 


OiCi 


Wi Wi 

1 a?6i7r^ 
^6 — 


- 2 




Wi 


wX 


--a>2-2 2 


2_^ 
3 wi 
02^2 


+2 


O2C2C2 2 af 1 a^b^T^ 


+ 


W2 

poo 


3 W2 


■ r d f 


U2UI 


UiU*2 


W2 

1 ay}2'K'^ 

6 Wn 


dx. 


PL 


a|c| 

W2 


Although the effective Lagrangian is expressed in 1,2 co¬ 
ordinates (for brevity), the Lagrangian must be expanded 
into the ± coordinates in order to evaluate the physical 
limit. For all the parameters we made the following ± 
coordinate substitutions into the expression for the total 
(effective) Lagrangian: 


Pi 


(2p+ +P-) 
2 


{2p+-p_) 

P2 = -;;- 


(29) 


c = 0, 
i = c, 
w = —2bw, 

0= -\w^ 


2^ • 


( 31 ) 


In the presence of the non-conservative term R, we expand 
in the ± coordinate systems and hnd the integrals: 


' —00 
poo 


' —00 
poo 


dn 

da_ 

dn 

WI 

dn 

dc- 

dn 


dx 


J PL 


dx 


PL 


dx 


PL 


d^- 

dn 

dw- 

dn 

d(t>- 


dx 


PL 


dx 


PL 


dx 


PL 


0 , 


TT^ea 


2 




0 , 

4ea^c 

w 


0 , 


4ea^ 

w 


Therefore, combining the conservative and non¬ 
conservative parts, the NCVA yields the following 
equations of motion for the NLS with linear loss: 


/ 

d = —at — ab, 
b = — 26^, 

TT^ TT^ ’ 

c = 0, 

^ = c, 
w = —2bw, 


= - W 


ic2 

O 1 


(32) 


with Pi e {ai,hi,Ci,di,uJi,^i} and P 2 G 
{02, 62, C2, (i2, W2j ^2}- 

Therefore, the full expansion with similar terms grouped 
together is Z = £1 — £2 + .R where 

d 

—— {u2ut — uiuX) dx. 

dp- J PL 

(30) 

Li and L 2 converge to the physical limit to recover the 
standard soliton evolution equations, £, i.e. the varia¬ 
tional approximation for the Hamiltonian, conservative, 
NLS equation. In this Hamiltonian case (i.e., in the ab¬ 
sence of perturbations), we obtain the following equations 



which correspond to the same dynamics as the conserva¬ 
tive case (1^ with the added loss term —at for the evolu¬ 
tion of the amplitude. 

The resulting dynamics comparison between the NCVA 
ODEs and the numerically integrated NLS is depicted in 
Fig. [T] for e = 0.1, (left set of panels) and e = 1 (right 
set of panels). The figure depicts in the respective top 
two panels the spatial profiles of the densities |up at the 
initial time (t = 0) and at a time of order 1/e for the 
PDE and ODE solutions and the remaining panels depict 
the evolution of the NCVA parameters. To compare the 
NCVA evolution of the parameters to the full NLS nu¬ 
merics, the numerical NLS solutions are projected (using 
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Figure 1: Evolution of an NLS bright soliton solution under the presence of linear loss of strength e = 0.1 (left) and e = 1 (right). A 
bright soliton, as described by Eq. isg, is used as an initial condition with the parameters: a(0) = u;(0) = 1, c(0) = 0.1, ^(0) = —5, 
and b(0) = 0(0) = 0. Depicted are the comparison of the NCVA approximation of Eq. II32II (red lines) with the full, numerical, NLS 
evolution of Eq. II21II (blue dots). The top two panels depict the density |up at the initial time (top subpanel) and at time t = 1/e (second 
subpanel). The bottom six subpanels depict the evolution of the NCVA ansatz parameters a, b, c, w, and tj>. For the full NLS evolution, 
the parameters are extracted by projecting the current solution into the NCVA ansatz using least-squares fitting. For e = 0.1 (left) the 
system is evolved for a total time of t = 1/e, while for e = 1 (right) the final time is taken to be t = 2/e. 


a least-squares fitting) onto the variational ansatz ua at 
discrete time intervals in order extract the variational dy¬ 
namical parameters associated with the entries of the vec¬ 
tor p. The time evolution for these parameters is depicted 
in the bottom six rows of panels in the figure. As observed 
from the hgure, the agreement between the NCVA system 
of ODEs and the numerical PDE integration is very good 
and accurately reflects the main dynamical features of the 
soliton solutions, such as the decrease of its amplitude (a), 
the increase of its width (hence the decrease of the inverse 
width controlling parameter w) etc. The high accuracy of 
the NCVA results, even in the case of large dissipation (see 
right set of panels in Fig. [T] where e = 1), may be deemed 
reasonable to expect in this case of linear dissipation. 

3.2. Density dependent loss 
As a second example, we use the attractive NLS equa¬ 
tion with a density dependent (nonlinear) loss term of 
strength e: 

iut + + \u\'^u = -i€\ufu. (33) 


to obtain the non-conservative ones. Taking R and ex¬ 
panding in the ± coordinates yields the non-conservative 
terms: 



on 


da- 

dH 

W 

dil 


J PL 


J PL 


dc- 

dn 


J PL 




PL 


dii 

dw- 

' dH 
d(f>- 


J PL 


PL 


dx 

dx 

dx 

dx 

dx 

dx 


0 , 

27r^e 4e 

9 3 

0 , 

8e a^c 
3 w 

0 , 

8e 

3 w ' 


We follow the same procedure as in the previous example. By combining conservative and non-conservative contribu- 
The conservative terms are the same and we just need tion, the NCVA yields the following equations of motion 
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Figure 2: Evolution of an NLS bright soliton solution under the 
presence of density dependent loss of strength e = 0.1. The full 
numerical solution is obtained from Eq. (I33I I (blue dots) while the 
NCVA results are obtained from Eq. (|34|l (red lines). Same layout 
of the panels as in previous figures. 


for the NLS with density dependent loss: 
a = — |ea^ — ab — 
h = 4«;4 _ 4 a 2^2 _ 2^,2 

TT-^ TT-^ ’ 


than an exponential). 


3.3. Exciton-polariton condensates 

The third, and hnal, example that we consider in the 
present work stems from the realm of exciton-polariton 
condensates. In this context, the condensing “entities” are 
excitons, namely bound electron-hole pairs. When con¬ 
fined in quantum wells placed in high-hnesse microcavities, 
these excitons develop strong coupling with light, forming 
exciton-photon mixed quasi-particles known as polaritons 
[ 3 ^ . An especially interesting feature of such polariton 
condensates is that their finite temperature leads the po¬ 
laritons to possess a finite lifetime —they can only exist 
for a few picoseconds in the cavity before they decay into 
photons. Hence, in this case, thermal equilibrium can 
never be achieved and the system produces a genuinely 
far-from-equilibrium condensate in which external pump¬ 
ing from a reservoir of excitons counters the loss of polari¬ 
tons. Exciton-polariton condensates offer numerous key 
features of the superfluid character for exciton-polariton 
condensates including: the flow without scattering (analog 
of the flow without frictio nl ISSll . the existence of vortices 
and their interactions [40ll4l| , the collective superfluid 


dynamics |42l | , as well as remarkable applications such as 
spin switches [^, and light emitting diodes M\ operating 
even near room temperatures. 

The pumping and damping mechanisms associated with 
polaritons enable the formulation of different types of 
models. One of these, proposed in Refs. 1^, sug¬ 

gests the use of a single NLS-type equation for the polari¬ 
ton condensate wavefunction which incorporates the above 
gain-loss mechanisms. Specifically, this model, based on a 
repulsive ((/ = —!) NLS equation with linear gain {ix{x)u) 
and density dependent loss {—i(T\u\^u) terrns, can be writ¬ 
ten in the following non-dimensional form 4^, 45| : 


c = 0, 

. (34) 

? = c, 

w = —2bw — %ea^w, 

y=1“^ “ 

which correspond to the same dynamics as the conser¬ 
vative case with the added nonlinear loss terms 

— (2/3 -I- 2/7r^)ea^ for the evolution of the amplitude and 
—d/TT^ea^ra for the evolution of the inverse width. The 
resulting dynamics for e = 0.1 are depicted in Fig. [2] We 
note that even larger values of the nonlinear loss (cf. e = 1) 
yield similar results. As for the linear loss case, we find 
very good agreement between the full NLS dynamics and 
the NCVA results despite the fact that the loss is in this 
case nonlinear. While the main properties of the solution 
dynamics (e.g. the decrease of the amplitude and increase 
of the width) persist, there are also nontrivial differences 
from the case of the previous subsection (e.g. in that the 
decay of the amplitude here follows a power law rather 


iut + - \u\^u - V(x)u = i [x(a;) - a\u\'^] u, (35) 

where a is the strength of the density dependent loss and 
we consider the localized, spatially dependent, gain 

X{x) = a exp > (36) 

induced by a laser pump of amplitude a and width /3. In 
the model, we use a general quadratic potential of the form 

V{x) = (37) 

To apply the NCVA, we take again the ui = ua{x, t;pi) 
and U 2 = UA{x,t','P 2 ) ansatze now defined as a Gaussian 
of the form 

UAix^t^p) =ae~2^e*(^“ + 4 ')^ ^•33^ 

where the ansatz parameter p) = {ai,Wi,bi,(j)i) for i = 1 
and 2 represent, respectively, the amplitude, width, chirp 
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Figure 3: Evolution of the ground state of Eq. ll35ll in the presence of a linear spatially dependent gain iIMt with a = 2 and /3 = 2, and 
density dependent loss of strength cr = 0.37, as well as a harmonic potential m of strength H = y/2. To craft initial conditions with 
amplitudes below (left set of panels) and above the equilibrium amplitudes we first computed the steady state of the NLS II35II which, after 
projection, using least-squares fitting, into the Gaussian ansatz ll38ll yields the following initial parameters: amplitude: a(0) = = 2.6431 

(equilibrium), width: u)(0) = 1.5583, chirp: b(0) = —0.1563, and phase: 4‘{0) = 0.2415. Then, we used, for the left set of panels, an initial 
condition with a(0) = 0.6608 = ae/4, i.e., four times smaller than the equilibrium solution; while for the the right set of panels we used 
an initial condition three times larger than the equilibrium solution, i.e., a(0) = 7.9292 = 3ae. Depicted are the comparison of the NCVA 
approximation of Eq. II42I I (red lines) with the full, numerical, NLS evolution of Eq. II35II (blue dots). The top two panels depict the density 
|up at the initial time (top subpanel) and at time f = 50 (second subpanel). The bottom four subpanels depict the evolution of the NCVA 
ansatz parameters a, 6, w, and (j>- For the full NLS evolution the parameters are extracted by projecting the current solution into the 
NCVA ansatz using least-squares fitting. 


and phase of the ansatz solution. It should be clarified 
here that our aim in this case (and in selecting this par¬ 
ticular ansatz) is to characterize the breathing motion of 
a ground state inside the trap, rather than to characterize 
the translational dynamics of the wavefunction (the latter 
would require a different ansatz, parametrizing the wave- 
function also by its center position). According to the 
NCVA method the Lagrangian is £ = £i — £2 +'^, where 
the conservative parts of the Lagrangian correspond to 


{uiulj. - 

+ 2 “i.xp + 2 



-bV(a:)|uip, 

(39) 

(■U2'«2,t - U2U2,t) 

+2 + 2 “2 



-bV(a:)|u2p, 

(40) 


and TZ has the same type of density dependent loss and a 
linear gain (equivalent to the negative of linear loss) shown 
in the previous test cases. The non-conservative terms are 


defined as follows: 


n = Pu*_+V*u-, (41) 

= —ix{x) {u 2 ul — U 1 U 2 ) + icr[|uip {u 2 u\ — u^tti) 

+ |u2p {U2U*1 — U 2 U 1 ) + U2U2ulul — UiUiU^ul], 


For all the parameters we made the substitutions of ± 
coordinates into the expression for the total Lagrangian 
and from the £1 and £2 parts we recover the conservative 
Euler-Lagrange equations for a Gaussian ansatz with four 
parameters. From the non-conservative term R, we ex¬ 
pand in the ± coordinate systems and hnd the integrals, 
which are combinations of the integrals for linear gain and 
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Finally, combining non-conservative and conservative 
terms, the NCVA yields the approximate equations of 
(breathing) motion for the exciton-polariton ground-state 
condensate of the form: 


(n = cra^w^ I ajSaw^ 

U — g cru 4 u , 2 + 2/32 2 ( u ; 2 + 2 / 32 ) 3/2 

3V2 , 2%/2a/3^a _ , 

2 u.2+2/32 + (u,2+2/32)3/2 




2 (u;2+2/32)3/2 


I 3V2aP^a‘^w , r,„..L 

+ u;^+2/3a 


: ^ 2 _ 1 
^ 8 “ 


(42) 


The comparison between the NCVA ODEs and the nu¬ 
merically integrated NLS is depicted in Fig. [3] for initial 
conditions below (left set of panels) and above (right set 
of panels) the equilibrium for the NLS. Below and above 
equilibrium refers, respectively, to initial solution ampli¬ 
tudes below and above those theoretically predicted to 
be at equilibrium. In the figure we use the coefficients 
a = 0.37, a = 2, /3 = 2, and ft = which ensures that 
the state with no excitation (i.e., without a dark soliton) 
is stable (see Ref. [i^). The top two panels show spatial 
profiles of the densities |up at the initial time {t = 0) and 
a final time of t = 50 for the NLS and NCVA solutions. 
The remaining panels depict the evolution of the NCVA 
parameters. As before, to compare the NCVA evolution 
of the parameters to the full NLS numerics, the numerical 
NLS solutions are projected (using least-squares fitting) 
onto the variational ansatz ua at discrete time intervals 
in order extract the parameters p. As observed from the 
figure, the agreement between the NCVA system of ODEs 
and the numerical PDE integration is good, although of 
lower quality in comparison to the rest of our considered 
cases. From Fig.[3]it is clear that both the original NLS dy¬ 
namics and its approximation using the NCVA are in very 
good qualitative agreement and good quantitative agree¬ 
ment. The lack of a better quantitative agreement stems 
from the fact that the solution to the original NLS prob¬ 
lem is only approximately a Gaussian —cf. configuration 
discrepancy between the converged states in the the sec¬ 
ond panel of the left set of panels in Fig. Actually, the 


solution is only close to a Gaussian for small atom num¬ 
ber, while upon increasing the atom number it approaches 
the so-called Thomas-Fermi (inverted parabola) profile H 
In all cases (below or above the stationary steady state), 
the dynamics of the NCVA and NLS converge (in an os¬ 
cillatory manner) to their respective stable solutions and 
do so consistently with respect to each other. 


4. Conclusions &: Future Challenges 

In this work we have extended the non-conservative vari¬ 
ational formulation recently proposed by Galley [l|. This 
was originally developed for classical mechanics, namely 
for systems with few degrees of freedom (and subsequently 
extended, including in the form of a variational approxi¬ 
mation, to nonlinear Klein-Gordon models in 0 ) , to the 
nonlinear Schrodinger equation (a complex partial differ¬ 
ential equation, i.e., infinite number of degrees of free¬ 
dom) . By using this non-conservative approach on a suit¬ 
ably chosen ansatz, it is possible to reduce the original 
infinite-dimensional dynamics to a system of ordinary dif¬ 
ferential equations on the ansatz parameters. We show 
that the resulting non-conservative variational method for 
the nonlinear Schrodinger equation is equivalent to the 
(linear) perturbative variational method. We provide sev¬ 
eral examples to test the validity of the non-conservative 
variational approach. In particular, we include examples 
with linear and density-dependent (nonlinear) loss. We 
also showcase the application of this method to exciton- 
polariton condensates that are inherently lossy and need 
a pumping term to balance losses. In all cases we see 
a very good qualitative and also, in principle, quantita¬ 
tive agreement (with the partial exception of the exciton- 
polariton condensate) between the original dynamics and 
statics for the nonlinear Schrodinger equation and its cor¬ 
responding reduced non-conservative variational counter¬ 
part. This agreement seems to be preserved even when 
the non-conservative terms are of the same order of mag¬ 
nitude as the other (conservative) terms. This is in con¬ 
trast with perturbation methods that intrinsically rely on 
the non-conserved terms being small when compared to 
the conserved ones. 

It would be interesting to apply in more detail the non¬ 
conservative variational methodology studied here further 
to exciton-polariton condensates —more specifically— 
that are, intrinsically, open systems of high interest and 
significant impact to ongoing experimental efforts. In par¬ 
ticular, such an approach would be valuable towards de¬ 
tecting the boundaries for stability inversion reported in 
Ref. [481 ]. It was noted in that work that, surprisingly, for 
some parameters values, the (originally) “excited” dark 


^In principle one could use a better suited ansatz like the q- 
Gaussian proposed in Ref. [3 at the expense of obtaining more 
complicated reduced NCVA ODEs. 
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soliton state (with one nodal point) becomes stable in fa¬ 
vor of the nodeless cloud state (the state without the dark 
soliton corresponding to the original ground state of the 
system in the absence on non-conservative terms) which in 
turn loses its stability. A systematic study of the dark soli¬ 
ton and its stability in such condensates can be found in 
Ref. [H^. Also valuable would be to extend this approach 
to the two-dimensional case where the equivalent stability 
inversion for vortices and rotating lattices has also been 
reported [45| . 

This work opens, more broadly, a few interesting av¬ 
enues for future explorations. In particular, it opens the 
possibility to apply a similar methodology to other models 
described by partial differential equations which contain 
non-conservative terms. The non-conservative variational 
methodology could prove very useful in cases where tra¬ 
ditional perturbative techniques, relying on the smallness 
of the non-conserved terms, fail due to the magnitude of 
the perturbations. One such model is the quintessential 
complex Ginzburg-Landau equation [l^ that models an 
extremely wide range of open systems including nonlinear 
waves, phase transitions, superconductors and superfluids, 
among others. 
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